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In quantum mechanics, predictions are made by way of calculating expectation values 
of observables, which take the form of Hermitian operators. It is far less common to 
exploit non-Hermitian operators to perform measurements. Here, we show that the 
expectation values of a particular set of non-Hermitian matrices, which we call column 
operators, directly yield the complex coefficients of a quantum state vector. We provide 
a definition of the state vector in terms of measurable quantities by decomposing the 
column operators into observables. The technique we propose renders very-large-scale 
quantum states significantly more accessible in the laboratory, as we demonstrate by 
experimentally characterising a 100 000-dimensional entangled state. This represents an 
improvement of two orders of magnitude with respect to previous characterisations of 
discrete entangled states. Furthermore, in numerical studies, we consider mixed quantum 
states and show that for purities greater that 0.81, we can reliably extract the most 
significant eigenvector of the density matrix with a probability greater than 99%. We 
anticipate that our method will prove to be a useful asset in the quest for understanding 
and manipulating large-scale quantum systems. 


One of the current challenges in the field of computing is harnessing the potential processing power 
provided by quantum devices that exploit entanglement. Experimental research aimed at overcoming 
this challenge is driven by the production, control and detection of larger and larger entangled quantum 
states ms- However, the task of characterising these entangled states quickly becomes intractable as 
the number of parameters that define a many-body system scales exponentially with the system size. 
To keep up with the ever-growing quantum state dimensionality, much effort is put into developing 
efficient characterisation methods (& 19 . 


Quantum state tomography is the process of retrieving the values that define a quantum system. 
The process typically involves two steps: i) gathering an informationally complete set of data and ii) 
finding the quantum state most consistent with the data set using post-measurement processing such 

Many efficient tomographic methods capital- 


as the maximum likelihood estimation algorithm 20 


ize on the first step by making simplifying assumptions about the state 11 T9 
number of measurements required to uniquely identify it 


thus reducing the 


In particular, tomography via compressed 
sensing allows one to efficiently reconstruct quantum states based on the fact that low-rank density 
matrices, i.e. quasi-pure states, are sparse in a particular basis 15-17,21 . Compared to assumption- 


free tomography, compressive sensing provides a square-root improvement on the required number 
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of measurements 110 


This improvement enabled the reconstruction of the density matrices of a 6- 

the largest phase-and-amplitude measurement 


qubit state 16 and a (17xl7)-dimensional state 17 


of an entangled state reported to date. Although compressed sensing does not make use of maximum 
likelihood estimation, it does require non-trivial post-measurement processing. 


Recently, Lundeen et al. reported on the direct measurement of a wavefunction using a method 
that, for the first time, required no involved post-measurement processing 19 . Their method is 


based on weak measurements, whereby one weakly couples a quantum system to a pointer state and 
subsequently performs a few standard strong measurements on the pointer state. The outcome of a 
weak measurement is known as the “weak value”, and in the conditions exposed in Ref. 19 the weak 


value is proportional to a given state vector coefficient. The method of Lundeed et al. can be used in 
combination with the assumption that the quantum state at hand is pure, providing the same square- 
root improvement as compressed sensing. Variations on the original scheme allow measurements of 
mixed states and increased detection efficiency 
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An important contribution of the work by Lundeen et al. was to link the state vector elements to 
the expectation value of weak measurements. We take a different approach, and point out that the 
enabling feature that allows access to the complex state vector is not weak measurement but the use 
of particular non-Hermitian operators. Although weak measurements provide a way to decompose 
these non-Hermitian operators, it is not the only suitable approach. Moreover, the introduction of 
weak values in the measurement procedure adds complexity to the experiment and the formalism that 
links weak values to measurement outcomes involves an approximation that breaks down in a variety 
of circumstances 


24-26 


In this paper, we propose an alternative approach to the direct measurement of quantum states that 
is exact in the case of pure states, proves to be reliable in the presence of noise, and is consistent with 
results obtained with well-established tomographic techniques. The key principle of our formalism is 
to decompose the particular non-Hermitian matrices that yield the complex state vector coefficients 
using only observables. Our method therefore only requires strong measurements, as in standard 
tomography, while maintaining the directness of weak-value-assisted tomography. The simplicity 
in both the experimental procedure and post-measurement processing renders our method ideally 
suited for the characterisation of large-scale systems, which can be high-dimensional, many-body 
or both. We begin by developing the theory on which our method is based and then demonstrate 
the potential of this scheme by experimentally retrieving the complex coefficients of a (341x341)- 
dimensional entangled state. 


Consider a quantum system in a d-dimensional Hilbert space, whose state vector 


d-1 

l») = £c,-|j> (i) 

i=o 

is expanded in the basis {|j)} and where Cj are unknown complex expansion coefficients. In order to 
retrieve these coefficients, we introduce the column operators Cj = \a)(j\, where |a) is an arbitrary 
reference vector. Each column operator has an expectation value 

(Cj) = <*| a) Cj , ( 2 ) 
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which is proportional to a complex state vector expansion coefficient. Since the value of (\&| a) is 
independent of j, we can express the state vector in terms of the column operators up to a phase 
factor: 


l*> 


7 e». 

3=0 


( 3 ) 


where v = |(T|a)| is a normalization constant. We can ignore the phase factor e*^ since it bears no 
physical significance. 


Most column operators Cj are not Hcrmitian matrices and are thus not observables. To overcome 
this apparent constraint, we recognize that any non-Hermitian matrix can be constructed from a 
complex-weighted sum of Hermitian matrices. Hence, the crucial step to our method is to construct 
the column operators in terms of measurable quantities: Cj = Yh q w jq^j q i where Wj q are complex 

weights and O n are observables. As a result, this allows us to retrieve any state vector element with 
a complex-weighted sum of measurement outcomes: 

c j ~ ~ w jq(Ojq)- ( 4 ) 

<? 

Equation [4] is an exact definition of the pure state vector that is provided in terms of measurable 
quantities. The above formalism readily applies to a general class of quantum states, including high¬ 
dimensional and many-body systems. 


As an example, consider the case of a qubit |T) = c o |0) + Ci|l) with |a) = |0) as the reference 
vector. The first column operator C 0 is Hermitian and given by the projector 10) (01. The second 
column operator C\ = 10) (11 is not Hermitian but can be constructed a number of ways. The Erst 
construction - which, as pointed out earlier, is a key part of the weak value formalism - is the complex- 
weighted sum of Pauli matrices: C\ = (a x + ia y )/2, a decomposition that requires two observables, 
each of which is made of two projectors or eigenvectors. A second decomposition requiring only three 
projectors is given by 

s > = El e ‘ 2 ' 8/3 K>M' (5) 

<7=0 d 

where |s 9 ) = (|0) +e* 47rg/ ; | 1))/a/2 are the states onto which the observables 0\ q project. In both cases, 
the qubit state vector is exactly given by |T) = ((Co)|0) + (Ci)|l))/(Co) 5 . 


To demonstrate the power and scalability of our scheme, we apply it to the measurement of a state 
entangled in greater than 100 000 dimensions. We provide a complete characterisation of the spatially 
entangled two-photon held produced through spontaneous parametric downconversion (SPDC). In 
general, SPDC can give rise to spatial and frequency correlations between two photons (4,27- 


36]. 

The purity of the spatial part of the full state can only be guaranteed if the two types of correlations 
are completely decoupled, which can be achieved in the collinear regime 


27 - see Supplementary 


Information section A for a theoretical estimation of our system purity. The consequences of applying 
our scheme to a quantum state with non-unit purity, which is always the case in the presence of noise, 
will be discussed below. 
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We express the spatial part of the entangled state in a discrete cylindrical basis of transverse spatial 
modes. The azimuthal part of the modes is given by e“0 where £ is an integer between — oo and oo 
and 0 is the azimuthal angle. This type of phase profile is known to carry £ units of orbital angular 
momentum (OAM). We decompose the radial part of the held with the recently introduced Walsh 
modes, labelled by the integer k ranging from 0 to oo 
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The Walsh modes all have the same 
Gaussian amplitude envelope, but different 7r-steps radial phase profiles. Combining the OAM modes 
with the Walsh modes yields a complete basis for coherent two-dimensional images. To perform the 
characterisation of the two-photon spatial held, we consider 31 OAM modes and 11 Walsh modes for 
each photon. The state vector thus takes the form 


15 


10 -15 10 




( 6 ) 


<1=-15 fci=0 ^2 = 15 k 2 =0 


Using the column-operator decomposition described in the Methods section, we sequentially measure 
all 116 281 coefficients which are shown in figure [iji and 1 3. The total Hilbert space dimension¬ 
ality of this measured state is more than two orders of magnitude larger than any previously reported 
amplitude-and-phase-characterised discrete entangled state 17 . As a simple verification of the accu¬ 


racy of our method, we calculate the probabilities associated with each joint mode via the Born rule, 

\„h,k 2 12 

IGi.fcil 

shown in hgure [l^, showing that we retrieve the correct magnitude of the amplitudes. 


, as shown in figure [lji. This result is consistent with the directly measured correlation matrix 


To rigorously assess the validity of the directly measured complex quantum state 10), i.e. both 
the amplitudes and the phases, we compare it to the results obtained through full tomography (i.e. 
assumption-free tomography). As full tomography cannot be performed on a (341x341)-dimensional 
entangled state in a reasonable time, we characterise a (5 x 5)-dimensional subset of the SPDC two- 
photon state. We perform the comparison in a basis of various OAM modes (£\ £ {1, —1,2, —2, 3}, 
£ 2 £ {1, —1, 2, —2, —3} ) and a fixed radial Walsh mode (ki = k 2 = 0). The total number of un¬ 
known parameters in the corresponding density matrix is equal to 624. After performing the direct 
measurement procedure in this basis, we record 8000 random projective measurements that we break 
into 8 sets of 1000. For each set, we recover a density matrix p cxp and calculate its purity and the 
fidelity with the directly measured state \ip); fidelity is defined as a/( 0|p ex p|0)- On average, the purity 
calculation yields (0.96 ± 0.02), and the fidelity gives (0.985 ± 0.004), where the uncertainties corre¬ 
spond to one standard deviation. After reconstruction of a density matrix, we fold that the average 
error between the measured count rates and the count rates predicted by the density matrix is 5.5%. 
This can be explained by shot noise, the pixelated nature of the SLM, and the finite aperture of the 
optical elements. While we expect unit purity, the 5% noise level accounts for the discrepancy with 
the measured value. 


The extremely high fidelity between the tomography results p exp and the directly measured state 
|0) indicates the validity of our approach for quantum state measurements applied to near pure states. 
To evaluate our method in the context of mixed states, we perform a series of numerical simulations 
where we vary the rank, purity, and dimension of an unknown state p sim , where no sources of noise 
are added to the simulated measurement outcomes. We apply our direct measurement procedure to 
these states and calculate the fidelity |(0|0 sim )|, where |0 S im) is the eigenvector of p s j m with the largest 
eigenvalue. For initial states p S i m with purity greater than 0.81, we measure a fidelity greater than 
0.99 in at least 99% of the cases. The dependency of this result on the dimensionality of the state is 
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negligible. This result indicates that our direct method is able to extract the primary eigenvector of 
a density matrix, even for a partially mixed state. Full details of this analysis and the density matrix 
reconstruction are presented in the Supplementary Information. 


Knowledge of the amplitude and phase of the state vector elements allows us to perform otherwise 
inaccessible calculations. As an example, we perform a calculation of the Schmidt decomposition 38 . 
This is equivalent to the singular value decomposition for the case of optical transfer matrices. The 
Schmidt decomposition yields a new joint basis in which the photons are perfectly correlated and 
where the joint modes have equal phases, as shown in figure [l]L When the Schmidt decomposition is 
applied to the entire state, we calculate a number of Schmidt modes equal to 142; this represents the 
effective number of independent joint modes contained within the state (the maximum for a (341x341)- 
dimensional state being 341). The Schmidt decomposed two-photon held is a good candidate for the 
violation of very-high-dimensional Bell inequalities 29 . Further details on the Schmidt decomposition 
can be found in the Supplementary Information. 


There are a number of approaches to reducing the necessary cost and effort for measuring large-scale 
quantum states. These include, but are not limited to, developing technologies for mode sorting 


42 


and arbitrary unitary transformations 43,44 , reducing the required number of measurement settings, 
and circumventing the requirement for reconstruction procedures. It is clear that there is significant 
interplay between each of these approaches. The theoretical implementation of an approach that 
combines the principles of our work with generalised measurements, such as POVMs (positive operator 
value measures), is considered in the Supplementary Information. The ability to use POVMs in the 
laboratory relies on the aforementioned technologies. Access to these types of technologies would 
reduce the overall number of measurement settings to uniquely recover a quantum state. However, 
such a system requires arbitrary unitary transformations for spatial states, which is in itself an active 
area of research 42-44 . Given the limitations of mode sorters for very large dimensions, and the 


practical nature of projective measurements, our scheme provides a simple and elegant method for 
the characterisation of large-scale quantum states. 


Our scheme allows direct access to the complex coefficients that define large-scale quantum states. 
The main result of our work is a novel method for retrieving a state vector coefficient with a complex- 
weighted sum of strong measurement outcomes. One challenge in reconstructing a quantum state from 
measurement outcomes lies in data processing; our scheme trades the difficulty of data processing for 
theoretical analysis prior to the experiment, that is, finding the measurements one has to perform. 
We anticipate that our work will have an impact on a number of disciplines, for example, quantum 
parameter estimation, measurement in quantum computing, quantum information and metrology. 


Methods 

Experiment The two-photon field is generated via SPDC with a 405-nm laser diode pumping a 
1-mm-long periodically-poled KTP (PPKTP) crystal with 50 mW of power. The experimental setup 
is shown in figure [2] We separate the two photons with a right angle prism and image the plane of the 
crystal to a Holoeye spatial light modulator (SLM) with a magnification of -10. We simultaneously 
display two holograms, one on each side of the SLM, to control the amplitude and phase profiles of the 
two photons independently. In order to make projective measurements of superposition modes, we 

We image the plane of the SLM with a magnification of —1/2500 


40 


make use of intensity masking 
to two single mode fibers. The combination of the SLM and singles mode fibers allows us to make 
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arbitrary projective measurements. All measurements are performed in coincidence with two single 
photon avalanche detectors, with a timing window of 25 ns, an integration time of 1 s for modes out¬ 
side the diagonal and 20 s for the diagonal elements (£i = —£2 and k\ = A^)- We start an automatic 
alignment procedure with the SLM every four hours to compensate for drift. Including the time it 
takes to calculate and display a hologram (about one second), the entire experiment takes two weeks; 
assumption-free tomography would take more than four centuries at the same acquisition rate. We 
perform no background subtraction and use the fundamental mode {£\ = —£2 — k\ — k 2 — 0) as the 
reference vector |a). The count rate of the fundamental mode is approximately 900 coincidences per 
second and varies by 10% over 24-hour periods. To correct for long term drift, we normalise each 
outcome to the count rate of the fundamental mode, which we measure before the measurement of 
each column operator. In standard tomography, the calculation of error bounds on the measured state 
is not a straightforward task 411. Here, we can calculate the error bound on a given coefficient with 


a weighted sum of the detector counts used to retrieve it. For a given state vector coefficient, the 
errors on the amplitude \cj\ and phase arg(cj) are both inversely proportional to the overlap v of the 
reference vector with the quantum state. In order to minimize the errors, it is important to choose a 
reference vector that has a high probability of occurrence within the state - the fundamental mode is 
the most probable one in our case. 


Two-body column-operator decomposition In order to decompose a given state vector coeffi¬ 
cient c/'fo into a set of measurement outcomes, we need to find a projector decomposition of the 


corresponding column operator = 10, 0)(£ 1 , Aq| ® 10, 0)(£2, k 2 1, as in equation 4 We numerically 


fold this column-operator decomposition, i.e. the complex weights w q and the observables O q , using 
the differential evolution algorithm (see Supplementary Information part D). By inspection, we find 
that the corresponding analytical form of the state vector coefficients is given by 


Pi(2,k2 


l £ jWVirfXoi* 

<7=0 


|S2, 9 )(S2, 


g\i 


( 7 ) 


where a/ 2|s m ,<j) = |0, 0)+e l47r9// '| £ m ,k m ) with m = {1,2}, and u = |(T|0,0)| is a normalisation constant. 
This decomposition is only valid when the state of any photon is different from the reference vector, 
i.e. \£m,k m ) % |0,0). Each coefficient measured with the above column-operator decomposition 
requires five projective measurements, thus explaining the 5 D 2 scaling, where D is the Hilbert space 
dimensionality of a single particle. The protocol scales much more favorably than assumption-free 
tomography, which requires D 4 projections. 


Here, we briefly explain our protocol for measuring the entire SPDC state vector. We measure 
more than 99% of the coefficients using the decomposition of equation [ 7 } The remaining column op¬ 
erators are the special cases 10, 0)(£i, Aq \ ® 10, 0)(0,01 and |0, 0)(0, 0| <g) 10, 0)(£2, ^ 2 1, which respectively 
correspond to a row and a column of the result shown in figure [l^. These column operators can be 
decomposed into only three joint local measurements using the projector 10, 0) (0,01 on one system 
and a column-operator decomposition similar to that of equation [5] on the other system. Finally, the 
column operator 10, 0) (0, 01 <g) |0, 0) (0, 0| is a projector, and its expectation value can be measured in 
a single experimental configuration. 

Full quantum tomography We perform full tomography with high count rates in order to achieve 
high accuracy. We set the magnification between the plane of the SLM and that of the single mode 
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fibers to 1/400. In this condition, we obtain a count rate of approximately 18,000 counts per second 
for the fundamental mode and integrate over 1 second for each individual projective measurement. 
The increase in the count rate of the fundamental mode comes at the price of lower count rates for 
high order modes. Regarding the full tomography measurements, we take an overcomplete set of 
1000 random projective measurements in a (5 x 5)-dimensional space. To minimize high-frequency 
components on the SLM, we limit the random superpositions to two-dimensional subsets of the state 
space. 
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a Amplitude (A) and phase of the measured state vector coefficients c e 2 ’^ 2 


Phase (rad) 




d Schmidt decomposition 


c Probabilities 


J-iM 


calculated from b 


e Directly measured probabilities 


b Renormalised (33 x 33)-dimensional subset 



Figure 1: Measured and calculated properties of the two-photon state, a We illustrate the complex 
state vector coefficients in matrix format. This representation is similar to that of an optical transfer matrix, 
where the lateral axes correspond to input and output modes. Here, each lateral axis corresponds to the 
spatial state of one photon of the pair. The OAM values l range from -15 to 15. For a fixed OAM value, 
the radial index k ranges from 0 to 10, thus the combined state has 341 x 341 dimensions. The amplitude 
and phase values of a coefficient are given by the height and color of a bar, respectively. For clarity, we 
darken the off-diagonal part. The small subset b of the state shows the phase gradient across the diagonal 
elements, which is typical of a Gouy phase shift, a common property of light passing through focus 45 . 
The corresponding calculated probability matrix c is consistent with the directly measured probabilities e. 
Finally, we calculate the Schmidt decomposition d of subset b, which gives the joint basis in which the subset 
can be expressed with the lowest number of modes. The indices S i and S 2 correspond to the states of each 
photon in this basis. 
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Holograms for the 

column-operator decomposition of C' 2 ^* 



Figure 2: Generation and characterisation of a two-photon field. The entangled state is pro¬ 
duced via SPDC in a PPKTP crystal and spatially separated by a prism. For the state determination 
stage, the crystal plane is imaged onto a spatial light modulator (SLM), which is in turn imaged to 
the input facet of two single mode fibers. In order to make a given projective measurement, we display 
the corresponding joint mode on the SLM and measure the coincidence rate between the two single 
photon avalanche diode detectors. The inset shows the five joint holograms that correspond to the 
column-operator decomposition of C' 2 o° . The state vector coefficient c 2 2 q° is given by ^ ^ (? ((9, ? )e* 47r9//5 , 
where the expectation value of a given observable is proportional to the measured count rate when 
displaying the corresponding hologram. 
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